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Abstract 

We argue the exact universal result for the three-point connectivity of critical percolation 
in two dimensions. Predictions for Potts clusters and for the scaling limit below p c are also 
given. 



1. The connectivity functions P n (x±, . . . ,x n ), i.e. the probabilities that n points belong 
to the same finite cluster, play a fundamental role in percolation theory [HE]- I n particular, 
their scaling limit determines the universal properties that clusters exhibit near the percolation 
threshold p c . Although most quantitative studies focus on the two-point connectivity, which 
determines observables as the mean cluster size S = P2( x ,0), also the connectivities with 
n > 2 carry essential information about the structure of the theory. In this paper we consider 
the case n = 3 for clusters in random percolation and in the g-state Potts model {q < 4) in two 
dimensions, in the scaling limit on the infinite plane. In particular, we will determine exactly 
the universal quantity 

o/ \ ^3(^1,^2,^3) m 

R(Xl,X2,X 3 ) = —= (1) 
y / P2(x 1 ,X 2 )P2(xi,X3)P 2 (x2,X3) 

at p c , as well as some of its features below p c . The connectivity combination (pQ) was studied in 
[3 H] for the case in which two of the three points are located on the boundary of the half-plane, 
it was shown to be a constant at p c and determined numerically and also analytically, exploiting 
the fact that points on the boundary yield linear differential equations for the connectivities 
at criticality [5]- It is not known how to write differential equations for random percolation 
connectivities if all the points are in bulk. 

It is well known that random percolation can be described as the limit q — >■ 1 of the g-state 
Potts model with lattice hamiltonian [6j [7] 

H = -J 6 s(x),s(y)> s(x) = l,...,q, (2) 

<x,y> 

as a consequence of the fact that the partition function Z = e ~ U can ^ e rewritten as [8] 

z = J2p nb (i-p) nb q N % (3) 

G 

where the sum is taken over all graphs G obtained putting bonds on the lattice, is the number 
of bonds in the graph, the number of absent bonds, p = 1 — e~ J the probability that a bond 
is present, and iV c the number of connected components (clusters) in G, with isolated sites also 
counted as clusters. Notice that ([3]) provides a continuation of the Potts model to real values 
of q. The random percolation transition corresponds to the limit q — > 1 of the ferromagnetic 
Potts transition, which is second order as long as q does not exceed a value q c (q c = 4 in two 
dimensions [9]). The connectivity functions can be related to the correlation functions of the 
Potts spin variable a a (x) = 5 s ( x ^ a — 1/q, a = 1, . . . , q; for p < p c , namely when the probability 
of finding an infinite cluster is zero, the relation for the two- and three-point functions reads 

q 2 

P 2 (x 1 ,X2) = -^—j (a a (xi)a a (x2)} , (4) 
q 3 

P 3 (X1,X2,X 3 ) = _ _ 2 ^ {(T a (xi)a a (X2)(Ta(x3)) , (5) 

where for random percolation the limit q — > 1 is understood in the r.h.s. It can be shown [10] . 
however, that these relations, evaluated for a generic q, give the connectivities for the clusters 
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entering the expansion ([3]), the so called Kasteleyn-Fortuin (KF) clusters (or "droplets" [I]) of 
the Potts model. 

It will be relevant for our purposes that in two dimensions the Potts model exhibits a duality 
property [6] relating a value K = e J — 1 of the coupling to the value K* = q/K. The self-dual 
point K c = y/q is the critical point. Duality allows to relate the correlation functions of the 
spin variables computed at a value J of the coupling to the correlation functions of "disorder" 
variables n a p, a, /3 = 1, . . . , q, a ^ f3, computed at the dual value J*. In particular one obtains 
the relations 

((T a (xi)a a (x 2 ))j<j c = -^(q - 1) (v a p{xi)np a (x2))j* , (6) 
(a a {x 1 )a a (x 2 )a a (x 3 )) J < Jc = -^(q - l)(q - 2) (/x a/3 (zi)/^ 7 (x 2 )/i 7a (z 3 )} j* . (7) 

Comparison with (Jl|5]) shows that P 2 and P 3 are related to the disorder correlators without any 
q-dependent prefactor. 

In the scaling limit {q < q c , J — > J c , distances \x{ — Xj\ = rij much larger than the lattice 
spacing) that we consider from now on, the lattice variables become fields of the Potts field 
theory, i.e. the simplest field theory realizing the S q permutational symmetry characteristic of 
the hamiltonian ([2]). Up to overall constants depending on the arbitrary normalization of the 
field a a (x), the correlators in ([MS]) are universal functions of the distances measured in units 
of the connectivity length. The normalization constants cancel in the combination ([1]), which is 
then completely universal. 

2. At criticality, conformal invariance implies that the two- and three-point correlators of 
scalar fields Ai(x) with scaling dimension Xi have the form [11] 

(A 1 (x 1 )A 2 (x 2 )) = S Xl , X2 C 12 , (8) 

(A 1 (x 1 )^ 2 (x 2 )A 3 (x 3 )) = C 123 r ^- Xl ~ X2 rfr Xl ~ X3 ^r X2 ~ Xi ■ (9) 



If the fields are chosen and normalized in such a way that Cy = 5ij, the constants Cjjfc, which are 

~<k 

'J 



invariant under permutations of their indices, coincide with the coefficients of the operator 



product expansion (OPE) 

Mxi)M*2) = C ij r^- Xl ~^ A k {x 2 ) + .... (10) 

Denoting X a the scaling dimension of the Potts spin field a a (x), we see recalling (j3H5j that 
P 2 oc r~ 2X<T , P 3 oc (ri 2 ri 3 r 23 )~ XtT , and 

q^ 

R(x 1 ,X2,x 3 ) = R c = _ 2 ^ C aaaaaa , p = p c , (11) 

where we have chosen 

C ~ ^ (12) 
in such a way that, in particular, P 2 is finite at q = 1. 
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In the two-dimensional case, to which we restrict from now on, the critical scaling Potts 
model with q < 4 is a conformal field theory [T2] with central charge 

c=1 -i<rnp (13) 
^- 2 ^wtt;- (14) 



with t related to q as [13] 
In conformal theories with c < 1, the formula 



[(t + l)m - tn] 2 - 1 

determines, for m and n positive integers, the scaling dimensions of scalar "degenerate" primary 
fields having the property that correlation functions containing them satisfy linear differential 
equations of order mn. For some discrete values of c, and in particular for any integer t > 2, there 
exist "minimal models" whose space of fields decomposes into a finite number of subspaces, each 
one originating from a degenerate primary [12]; the number of primaries with a given scaling 
dimension is an integer in minimal models 1 15j . In the Potts model, the fields a a have |16U13| 

X a = X( t+ i)/ 2) (t+l)/2 (16) 

and multiplicity q — 1 (X) Q a a = 0), allowing minimality only for q = 2 (Ising model, t = 3) and 
q = 3 (t = 5); the leading S^-invariant field e, with X £ = X^^i, is degenerate for generic values 
of g. 

In [IT] the properties of the degenerate fields were exploited to compute the OPE coefficients 
Cf, for the "diagonal" series of minimal models, in which all the primaries appear with multi- 
plicity one. The Potts model belongs to this class for q = 2, but even for this case the results for 
the minimal OPE are not sufficient to determine R c . Indeed, it is clear from (|lip that a finite 
R c requires C (Ta(Ta(Ta = at q = 2, and this is trivially ensured by the spin reversal symmetry of 
the Ising model. Hence, even the determination of R c at q = 2 requires the computation of the 
OPE coefficient for continuous values of q\ even forgetting the problem with the multiplicity of 
the fields, the formulae of [17] can be evaluated only for the discrete values of t corresponding 
to minimal models. 

Few years ago Al. Zamolodchikov |I8j approached the problem of the derivation of the OPE 
coefficients of c < 1 minimal models within a conformal bootstrap method based on the use 
of correlators of four scalar fields, one of which being the degenerate primary with dimension 
X\ t 2 or X2,i ■ The other three fields are simply required to appear with multiplicity one; the 
mathematical treatment does not put any constraint on their scaling dimensions Xi, i = 1,2,3. 
The method leads to functional equations for the OPE coefficients which are related by analytic 
continuation to the functional equations arising in Liouville theory (c > 25) [19]. The solution 
for c < 1, however, is not an analytic continuation of the Liouville solution, and reads |18j 

C xl,x 2 = C x u x 2 ,x 3 = (17) 
A T(ai + a 2 -a 3 + /3)T(a 2 + a 3 - ai + /3)T(a 3 + ai - a 2 + /3)T(2/3 - + a x + a 2 + a 3 ) 



[T(2ai + /9)T(2 01 + 2/3 - /3" 1 )T(2a 2 + /3)T(2a 2 + 2/3 - /3" 1 )T(2a 3 + /3)T(2a 3 + 2/3 - /3" 1 )] 



i j 
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with 



0=y/t/(t+l), 

X i = 2a i {a i + P-r 1 ), 



.4 



T(/3) 



7(z) 



r(x) 
r(i-z) 



(18) 
(19) 

(20) 



T(x) = exp 




Q 



2 sinh^ 

t 



x e 



2 x ; 2 



sinh y sinh X 



Q = /3 + r 1 - (21) 



The integral in (|2ip is convergent for < x < Q; outside this range T(x) can be computed using 
the relations 



T(x + l//3) 



j(x/(3)/3 2x ^~ 1 T(x) . 



(22) 
(23) 



The function ()1T|) can be evaluated for continuous values of t and of the Xj's. Taken literally, 
it would give the OPE coefficients for a theory with arbitrary c < 1 and with a continuous 
spectrum of fields having multiplicity one. It is remarked in |18] that the consistency of such a 
theory, in particular from the point of view of modular invariance |144 115). is an open question. 
Concerning the values of c corresponding to minimal models, (|19p gives the scaling dimensions 
(fT5|) for Gtj equal to 

(n ~ W _ (m - l)r 1 (0A s 

Although a general proof is not available, checks case by case show that (I17j) reproduces the 
OPE coefficients of minimal models obtained in [17] . at least when these differ from zero. In 
some cases for which the minimal OPE prescribes vanishing coefficients, (|17j) gives instead finite 
numbers whose interpretation is considered "mysterious" in [18]. As an example, (|17p evaluated 
for X\ = X2 = X% = X2 } 2 at t = 3 does not vanish, despite the fact that the spin three-point 
function is zero in the Ising model. 

Recall now that the two-dimensional Potts model contains also disorder fields, dual to the 
spin fields and with the same scaling dimension X a . They satisfy the OPE 



Pap Hp-f = ?>aj(I + C e e + ...) + (1 — <5 Q:7 )(C At jU a7 + ...), 



(25) 



where we omit the coordinate dependence for simplicity and the coefficient in front of the identity 
is fixed to 1 by © and f)12|) . Use of ([6][7]) at the self-dual point then leads to the result R c = Cu- 
Sq-invariance gives to (|25[) a two-channel structure (a = 7 or a ^ 7) equivalent to that produced 
by two fields \x and jl satisfying 



I + C e £ + . . . , 

+ /!) + .... 



(26) 
(27) 
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The field (f>= (fi + fi)/ \/2 then satisfies 

<t>4> = I + C e e + -±L4> + ... , (28) 
V2 

namely a "neutral" OPE for fields with multiplicity one, as the one assumed in the derivation 
of (|17|) . Notice also that, due to the neutrality of the fields in (128p . has no reason to vanish 
for any value of q; in particular, this is not in conflict with the symmetries of the Ising model, 
because the absence for q = 2 of the vertex with three disorder fields in (|25[) is enforced by the 
factor 1 — 5 ai . 

These observations may suggest the following interpretation for the function (11T|> : it encodes 
the "dynamical" information about the OPE of conformal field theory for c < 1, and knows 
nothing about internal symmetries, which, on the other hand, are not uniquely determined by 
the value of c; symmetry considerations have to be developed separately and produce a dressing 
of (|17p by factors which may vanish, suppressing the vertices not compatible with the given 
symmetry. Letting aside the general validity of this interpretation, the OPE's (|25p and (|28|) 
make it plausible for the case we are discussing, and lead us to take 

Rc = V2C Xa>x „ >x „, (29) 

with X a given in (|16j) . The results for the integer values of q are given in Table [TJ the value 
for q = 4 is obtained in the limit t — > oo. In [20] the constant value ~ 1.022 was obtained 
numerically for critical percolation on a cylinder and ([T]) evaluated for x\ on one edge, X3 on the 
other edge and %i far away from both edges. Since in such a configuration the cylinder, seen 
from X2, looks infinitely long and becomes conformally equivalent to the plane, we see in this 
result a confirmation of our analytic value for R c at q = 1. 



q 


1 


2 


3 


4 


R c 


1.0220.. 


1.0524.. 


1.0923.. 


1.1892.. 


Rc 




1.3767.. 


1.3107.. 


1.1892.. 


1 KK 


1.0450.. 


1.1547.. 


1.3160.. 


1.8612.. 



Table 1: R c and are the values of the quantities (jlip and (|32p for KF clusters; q = 1 

corresponds to random percolation. R c is obtained from (I29p with X a replaced by X„, the 
scaling dimension associated to spin clusters. 

The result (|29p refers to KF clusters. Concerning the ordinary spin clusters, i.e. those 
obtained connecting nearest neighbors with the same value of the spin, they are also critical at 
J c in two dimensions [10], with connectivities related to the correlation functions of the field 
with scaling dimension = X t / 2t t/2 [ZD E2]. For the Ising case, the best understood field 
theoretically [23l [23] , the ratio (P) is expected to be given by (f29|) with X a replaced by X&. We 
give this value R c in Table [1] those obtained in the same way at q = 3, 4 are also quoted, but 
these cases are less clear. 



5 




Figure 1: Triangle identified by the points x±,X2,xs. The lines joining the vertices of the triangle 
to the Fermat (or Steiner) point Y form 2tt/2> angles. 



Comparison between (fTT]) and ([29)) determines C UaCTaCTa and its relation with C^. The re- 
lations C CTceCT/3 = (q5 a p - l)/q 2 and C CT1(TlCTl = (1 - q)C aiCTia2 = \{q - l)(g - 2)C (JlCT2<J3 are a 
consequence of the symmetry. Eqs. ([6]) and (f28l) lead to C £ = —Q 2 (q — ^) C% a a a = Cx a ,x a ,x e -, 
where we also took into account that e changes sign under dualitjo. The OPE between a a and 
/ig 7 produces parafermionic fields. An analysis similar to that of [25] suggests that they have 
spin this is certainly the case for q = 2, 3 |26| . 

3. Away from criticality the Potts field theory is solved exactly in the framework of the 
factorized S-matrix [27], from which large distance expansions can be obtained for the correlation 
functions |28[ [29] . The S q symmetry is more transparently implemented working in the ordered 
phase (J > J c ), where the elementary excitations (kinks) are interpolated by the disorder fields; 
the results for the disordered phase are obtained by duality. If xi,X2,x 3 are the vertices of a 
triangle whose internal angles are all smaller than 2ir/3, the asymptotic behavior of the correlator 
flTJ when all the distances between the vertices become large reads [27J ([30] for a derivation) 

Mxi)^(x2)^ a (x 3 ))j >Jc = -± T« K K (r Y /£) + 0(e-'/*) . (30) 

Here is the one-kink form factor of the disorder field, Kq(x) = dy e~ xcosh ^ y is a Bessell 
function, and ry = r\ + r2 + ^3 is the sum of the distances of the vertices of the triangle from 
the Fermat (or Steiner) point Y, which has the property of minimizing such a sum (Fig. [1]); 
ry < p = min{ri2 + 7*13, T\2 + ^23i r i3 + ^23}- £ is the connectivity length (inverse of the kink 
mass), also determined by the large distance decay of the two-point function 

Mxi)^ a (x2))j>j c = S K (r 12 /0 + 0(e~ 2r ^) , (31) 
1 For example ITT]) gives C e | q= 2 = -1/2. 
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and is the three-kink vertex given by [2 



1 KK 



1 . 2vrA 

- sin — - flf(AJ , 



g{\) = exp 



dt sinh I - 
3 



A sinh _ sinh 



t sinh cosh | 



(32) 
(33) 



with q = 2sin^j-, A G (0,3/2). It follows from ([M7]) that (I30ti3ip are precisely the asymptotics 
for large separations of the connectivties (jHE]) for KF clusters below p c . It follows 

K (r Y /0 



R(xi,x 2 ,x 3 ) ~ v^T^ 



Pc 



(34) 



The dynamical information is entirely contained in the three-kink vertex, whose values for q 
integer are giver§| in Table [TJ In the opposite limit, in which all the distances are much 
smaller than £ (always remaining much larger than the lattice spacing), the function R tends to 
its constant critical value R c . 

Finally, consider the ordinary spin clusters for the Ising model. It was argued in [23] that in 
this case the scaling limit for p —> p~ (in zero magnetic field) corresponds to a renormalization 
group trajectory with infinite connectivity length ending into a random percolation fixed point at 
large distances. Then one expects ([TJ to interpolate from R c \q=2, when all the distances between 
the points are small, to R c \q=i, when all of them are large. Although universal, this crossover 
could be easier to observe on the triangular lattice, which has p c = 1/2 for site percolation 
and is expected to minimize the corrections to scaling [23]. If instead we work at the Ising 
critical temperature and add a magnetic field H, the connectivity length is finite and the large 
separation behavior of (pQ) for clusters of positive spins and H — > 0~ has again the form (]34p . 
with T^ K = 5.7675.. given by (|32~j) evaluated at A = 5/2 [24]. Since the integral in (j33j) diverges 
for A > 3/2, one uses the analytic continuation 



5(A) 



r(i-f )r(i + A) 
r(i + |) 



exp 



00 ft 
dt sinh ( — 



sinh[| (l-i)]- e -* sinh + 



t sinh ^ cosh | 



for 3/2 < A < 3. 
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